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Abstract. This article summarizes our understanding of the Kondo effect in 
graphene, primarily from a theoretical perspective. Wc shall describe different ways 
to create magnetic moments in graphene, either by adatom deposition or via defects. 
For dilute moments, the theoretical description is in terms of effective Anderson or 
Kondo impurity models coupled to graphene's Dirac electrons. Wc shall discuss in 
detail the physics of these models, including their quantum phase transitions and the 
effect of carrier doping, and confront this with existing experimental data. Finally, 
we point out connections to other quantum impurity problems, e.g., in unconventional 
superconductors, topological insulators, and quantum spin liquids. 
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1. Introduction 

The low-temperature behaviour of dilute magnetic impurities in metals, known as 
the Kondo effect, is a prime example of electron-correlation physics. The impurity's 
magnetic moment is screened below a temperature Tk by the formation of a many-body 
singlet with the conduction-electron bath. The Kondo temperature Tk itself depends in 
a non-analytic fashion on the Kondo coupling and the bath density of states, signalling 
the breakdown of perturbation theory. Starting with Kondo's work in the 1960s [1], 
we now have an essentially complete set of theoretical descriptions of the Kondo effect 
available [2], and agreement between theory and experiment has been established on a 
quantitative level. 

The understanding of the Kondo effect in metals has prompted to investigate similar 
impurity physics in other settings, with the overarching goal to employ impurities as 
local probes of the host's properties . This article is devoted to a particularly interesting 
and timely case, namely the Kondo effect in graphene. Here, the impurity spin interacts 
with the Dirac fermions of the two-dimensional (2d) sheet of carbon atoms [3, 4, 5, 6]. 
For a local magnetic impurity in charge-neutral graphene, this results in the Kondo 
effect being qualitatively different from that in conventional metals, because the bath 
density of states now vanishes at the Fermi level: Kondo screening is suppressed at small 
Kondo couplings, and a non-trivial impurity quantum phase transition [7] between an 
unscreened and a screened impurity spin obtains. In fact, this phase transition has 
first been discussed for magnetic impurities in unconventional superconductors, where 
field theories and numerical solutions for the resulting pseudogap Kondo problem have 
been worked out. Graphene offers the attractive possibility of tuning the chemical 
potential relative to the Dirac point, such that the crossover between pseudogap and 
conventional Kondo physics can be accessed in detail. Furthermore, the 2d nature of 
graphene naturally allows one to employ scanning-probe techniques to locally study 
impurity physics. 

In this article we review the theoretical understanding of the Kondo effect in 
graphene, together with the current status of experiments. As we will discuss, clear- 
cut experimental verifications of some of the exciting theoretical ideas are lacking, and 
we contemplate on possible sources of comphcations. We also highlight theoretical 
connections between the Kondo effect in graphene and other quantum impurity 
problems, such as impurities on the surface of topological insulators. 

1.1. Outline 

The body of this article is organized as follows: In Sec. 2 we introduce the electronic 
structure of graphene and review general aspects of magnetic moment formation. 
We then discuss various possibilities of experimentally realizing magnetic moments 
coupled to graphene sheets, together with the relevant microscopic descriptions. Quite 
generically, this will lead to versions of the peudogap Kondo model, whose theoretical 
treatment is discussed in some detail in Sec. 3. We describe the phase diagram. 



Kondo Impurities in Graphene 



3 



the quantum field theories and resulting critical properties, the crossovers for finite 
chemical potential, and the implications for Kondo physics in graphene. Sec. 4 confronts 
these theoretical results with experimental data obtained on impurity-doped graphene. 
Various real-world complications and their infiuence on the interpretation of experiments 
will be discussed as well. Finally, in Sec. 5 we discuss quantum impurity problems which 
are relatives of the graphene Kondo problem, thereby highlighting the generality of the 
theoretical concepts developed in the field. 



2. Magnetic impurities in graphene 

In this section, we discuss how to realize magnetic impurity moments coupled to 
graphene conduction electrons. To set the stage we first summarize basic aspects of 
both the electronic structure of graphene and the formation of local moments in general. 



2.1. Electronic structure of graphene 

Graphene is a 2d hexagonal arrangement of carbon atoms [3, 4, 5, 6]. While the sp^ 
orbitals hybridize to yield the a orbitals which are electrically inert and responsible for 
the remarkable mechanical robustness of graphene, its electronic structure is determined 
by the orbitals which form the vr-bonds. This allows electron hopping between 
adjacent atoms and gives rise to a kinetic energy described by the following tight- binding 
Hamiltonian [8]: 

^0 = E {alA, + h.c.) - t' {al,a^^^ + b^X, + h-c) (1) 

Here, a^j^i and b^j are annihilation operators for electrons on sites i,j of the two 
sublattices A and B, see Fig. 1 a. The sums run over pairs of nearest and next-nearest 
neighbours, respectively, with the hopping matrix elements given hy t — 2.8 eV and 
t' ~ O.lt. Diagonalizing the Hamiltonian (1) yields two dispersive bands 



Ei,k = ± ty/SW^ - t'f^. With 

fj^ = 2 COS (ySkyUi^ 4cos ^"^^2/'^i^ cos ^k^ai^ , (2) 

where ai = 1.42 Adenotes the bond length, i.e. the distance between neighbouring 
carbon atoms. The two bands (dubbed n* and tt) touch at the two inequivalent 
wavevectors K = (g, ^) and K' = (g,-^), see Fig. la. Close to K and 
K' the dispersion is found to be hnear: 

E^,,^=±vf\^ (3) 

where q = k — K {oi K'), vp = 3tai/2 a; 1.1 x lO^m/s, and an additional constant 
proportional to t' has been omitted in E±. 

For charge-neutral graphene, the electronic system is half-filled, i.e. in the ground 
state all E_ states are filled while the £"+ states are empty, such that the Fermi level 
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Figure 1. a.) Honeycomb lattice of grapiiene with two inequivalent carbon atoms per 
unit cell, A and B, and its hexagonal Brillouin zone. The points K and K' are the 
touching points of the E± bands, b.) Sketch of the AIM: A local spin-degenerate level 
with energy is hybridized via v with a sea of conduction electrons. Local double 
occupancy costs the Coulomb energy U. 



coincides with the band energy at K and K'. Then Eq. (3) is the dispersion of the 
system's low-energy excitations which admit a description in terms of two-component 
massless Dirac equations, one for each of the valleys at K and K' [5]. The 2d linear 
dispersion results in a low-energy density of states (DOS) per spin which vanishes linearly 
at the Fermi level, 

2 , , 



\UJ\ 



(4) 



rendering graphene a semimetal. (Here and in the following, energies are measured 
relative to the Fermi level, unless otherwise noted.) 



2.2. Local-moment formation in metals 

The general mechanism of local-moment formation in metals was formulated by 
Anderson in 1961 [9]: a strong Coulomb interaction U between electrons in a spin- 
degenerate doublet of levels, can freeze out charge fluctuations, leaving behind an 
effective spin degree of freedom interacting with the spin density of the conduction 
electrons, . Typically, such a situation is realized for impurity atoms with partially 
filled dor f shells. The corresponding minimal model is known as the Anderson impurity 
model (AIM): 

where Uda = d'^^d^. Moment formation can be understood starting from the atomic limit, 
= 0. For ea < Ep and + f/ > Ep the d level prefers single occupancy, such that 
charge fluctuations are frozen out and an effective spin 1/2 degrees of freedom remains. 
Upon switching on the hybridization f the so-formed local moment becomes entangled 
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with the conduction electrons. It is convenient to convert the momentum dependence 



which fully characterizes the impurity's bath. For small v^, one can utilize a Schrieffer- 
Wolff transformation to derive an effective Kondo model from Eq. (5), describing the 
interaction of the local-moment spin with the conduction electrons, see Sec. 3. 

2.3. Magnetic adatoms on graphene 

We now turn to the graphene-specific discussion of how to realize magnetic impurity 
moments with sizeable electronic coupling to the host electrons. An apparently 
straightforward route is to place a magnetic ad-atom onto the graphene sheet, e.g., 
using the manipulation capabilities of a scanning tunneling microscope (STM). For 
magnetic atoms like Fc or Co on the surface of conventional metals, this route has been 
successfully used in the past to study local spectral signatures of Kondo screening using 
STM techniques [10]. 

In the following we discuss theoretical aspects of such adatoms on graphene; 
experiments will be reviewed in Sec. 4. The key questions for a quantitative 
understanding of the adatom's magnetism are: 

(i) At which lattice position does the adatom adsorb? For graphene, possible high- 
symmetry locations are shown in Fig. 2a and labeled h (hollow, in the center of a 
hexagon), b (bridge, on a bond between two C atom), and t (on top of a C atom). 

(ii) What is the spin state of the impurity adatom? 

(iii) In an Anderson-model description, how are the impurity levels hybridized with 
conduction electrons? 

All these question turn out vital for the presence or absence of the Kondo effect. 
In general, answering these questions is highly non-trivial and requires ab-initio 
calculations, typically using variants of density-functional theory (DFT).J 

2.3.1. Co. For Co atoms on an isolated graphene sheet, DFT calculations using the 
generalized gradient approximation (GGA) found that the preferred adsorption is at site 
h with spin S — 1/2 [12]. However, upon accounting for the local Coulomb repulsion 
U and Hund's-rule coupling J within the GGA-I-U method, this picture was modified: 
While for values of C/ = 2eV and J = 0.9 eV the h position and S — 1/2 were still 
favoured, increasing U to 4eV selected the t position and S — 3/2 [12, 13]. As U can 
only be estimated to be in the range 2eV< U < 4eV, a clear-cut answer was missing 
here, but the case of a spin 5 = 1/2 on the h site was advertised as the most promising 
candidate for observing Kondo physics. 

I Notably, even for the "classic" situation of Fc atoms in gold or silver, the correct model description 
has only been determined very recently [11] to be a spin-| three-channel Kondo model. 



of and eg into an energy-dependent hybridization function 




(6) 
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In an alternative calculation, based on DFT augmented by a dynamical treatment of 
the 3d levels in the framework of the one-crossing approximation, dubbed GGA+OCA, 
a S = 3/2 configuration on the t position was found to be most stable [14]. Here, the 
authors argued in favour of a Kondo effect with full screening, as three conduction- 
band channels coupled to the impurity spin. Finally, in a refined quantum-chemical 
calculation [15] based on a complete active-space self-consistent field approach, albeit 
on small clusters, it was found that Co in an h position favours a higher-spin state of 
S^3/2. 

Overall, the situation concerning Co adatoms is unclear at present, and more theory 
work is called for. Specifically, the effect of the substrate, which possibly influences the 
adatom's behaviour, has not been investigated so far. 

2.3.2. NiH. Recently, it has been proposed that a more promising route towards 
realizing the Kondo effect could be provided by using NiH as adsorbing molecule [13]. 
From GGA-I-U, the molecule favors b. S = 1/2 state in the h position [13]. 

2.3.3. Effective model. As in Sec. 2.2, the physics of a locahzed impurity level 
hybridized with graphene electrons can be described by an Anderson impurity model, § 
which then may be mapped onto an effective Kondo model, see Sec. 3. 

For a spin- 1 Co atom in the h position, it was argued [12] that, due to an 

approximate orbital degeneracy, the impurity behaviour at elevated energies corresponds 
to that of an SU(4) Kondo effect. The orbital splitting is roughly 60meV, such that a 
standard single-channel SU(2) S* = 1/2 Kondo or Anderson model applies at energies 
below this scale, albeit with a non-standard hybridization function The latter 

was calculated in Ref. [12] and is reproduced in Fig. 2b. It vanishes near the Dirac-point 
energy (set to zero here) according to 

lmA{uj) oc |a;|, (7) 

a behaviour inherited from the graphene DOS, while for higher energies there is sizeable 
particle-hole asymmetry. Both features turn out vital for the Kondo effect, see Sec. 3.2. 
It can be expected on symmetry grounds that this effective model with a similar 
hybridization function applies to NiH on the h site as well. 

2.4- Defect-induced moments in graphene 

A different route towards magnetic impurities in graphene is via point defects which 
themselves induce moments. Relevant defects include vacancies, created e.g. by 
irradiation [18, 19], and hydrogen and fluorine adatoms. 

§ The mean-field solution of the Anderson model for charge-neutral graphene has been discussed in 
Ref. [16]; for aspects of the full solution see Ref. [17] and Sec. 3. 
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Figure 2. a) Possible high-symmetry adsorption sites of ad-atoms on the graphene 
sheet, denoted by h (hollow), t (top), and b (bridge), b) Hybridization function of 
the active impurity level of a Co atom in the h position in a spin state of = 1/2 
[12]. Sohd: DFT result, dashed: tight-binding fit. El, E2 refer to different orbital 
configurations, with El being lower in energy (taken from Ref. [12]). 



2. 4-1- Tx-orhital magnetic moment. Removing a single site from the 7r-electron tight- 
binding Hamiltonian in Eq. (1) induces a single localized state at the Dirac-point energy 
for t' = 0, which becomes a quasi-localized resonance near the Dirac-point energy for 
non-zero t' [5, 20]. Consequently, it has been proposed that, upon including Coulomb 
interactions, a magnetic moment may be formed in this localized state, i.e. in the 
vicinity of the vacancy. However, in real graphene the lattice near the vacancy will 
reconstruct, such that determining the proper effective model requires ab-initio studies. 

Notably, the results of those studies are again controversial. Initially, local-moment 
formation for a vacancy was confirmed using DFT, however, the coupling to the 
conduction electrons was argued to be ferromagnetic [21], such that no Kondo effect 
can be expected. A similar conclusion was drawn from a study using dynamical mean- 
field theory (DMFT) [22] where a Curie-type susceptibility was found indicative of a free 
moment. In contrast, for an H atom adsorbed on top of a C atom, Ref. [23] argued that 
the physics can be described by a single S* = 1/2 impurity coupled antiferromagnetically 
to the environment, such that screening with sizable Kondo temperature should be 
possible. Finally, a more recent DFT calculation [24] for a finite concentration of 
H adatoms concluded that, while there can be localized cr-orbital moments, 7r-orbital 
moments only occur for the unlikely situation of hydrogenation of all dangling a bonds. 
This paper pointed out the importance of considering both a and tt bonds and their 
reconstruction around a vacancy. 

In any case, if a vr moment forms then it can be expected that its coupling to 
the conduction electrons is described by an Anderson/Kondo model with a pseudogap 
hybridization function as in Eq. (7), possibly with a large particle-hole asymmetry at 
higher energies due to potential scattering. The physics of this pseudogap Kondo model 
will be discussed in detail in Sec. 3. 
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2.4-2. a -orbital magnetic moment. An interesting alternative is to consider the 
carbon's a orbitals. While local-moment formation driven by Coulomb interaction is 
possible here as well, the obstacle is that - for flat graphene - the hybridization between 
a and vr orbitals vanishes, i.e. such a moment would not couple to the conduction 
electrons. However, hybridization of a and tt orbitals becomes possible once structural 
deviations from the flat geometry are included, i.e. by corrugations of the graphene 
sheet around the impurity site. 

Ref. [25] investigated vacancies and cr-orbital moments from an LDA+U 
perspective. It was found that S — 1/2 or S — 1 a moments emerge, which can have 
sizeable single-channel hybridization to the tt electrons upon including rippling, which 
occurs under small isotropic compression around reconstructed vacancies. Interestingly, 
the hybridization function within an effective Anderson model was found to be low- 
energy divergent according to 

ImAHoc I , (8) 

\u\m \uj/D\ 

where D is the bandwidth, providing a high-energy cutoff. In that situation the Kondo 
temperature can largely be enhanced due to the massive density of states at low energies 
[25, 26], but detailed studies of this model are not available. 



3. The pseudogap Kondo problem 

We now discuss the rich physics of the so-called pseudogap Kondo model, relevant to 
low-energy behaviour of magnetic moments in graphene. We will restrict our attention 
to the case of a spin S = 1/2 coupled to a single screening channel; the two-channel 
version will be briefly mentioned in Sec. 4.2. 



3.1. The pseudogap Kondo model 

The standard Kondo Hamiltonian [2] reads 

^ = E 'A.'k. + ^0 E clc,,^ + JoS- .-0, (9) 

k,<T k,k',<T 

where the notation follows Sec. 2.2, S is the impurity spin S = 1/2, and sq = 
2 Sfcfc' ^\ '^'^^'^k' ' conduction-electron spin density at the impurity site, with 

r the vector of Pauli matrices. The Kondo coupling Jq and the potential-scattering 
strength Vq characterize the impurity. If the Kondo model is derived from the more 
general Anderson model, Eq. (5), in the limit of small charge fluctuations, second-order 
perturbation theory yields [2]: 

\\€d\ \U + €d\^ 2 \\€d\ \U + €d\^ 

where v^ = v has been assumed. Note that T^o 7^ breaks particle-hole symmetry if the 
host DOS is particle-hole symmetric, = p{—uj). 
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For a metallic host, the DOS p{uj) is finite at the Fermi level. Then, for 
antiferromagnetic J > the impurity spin is screened below the so-called Kondo 
temperature Tk. For a flat conduction-band DOS, p{uj) — po, one flnds [2]: 

Tk^ ^/DJoe-^/^'^°''°K (11) 

Importantly, in this metallic Kondo problem, the crossovers at finite energies and 
temperatures are characterized by the single scale 7k- For instance, the impurity 
susceptibility displays single-parameter scaling: Ximp(T) is a universal function of T/Tk 
only, and does not depend on further microscopic details. 

In the following, we will instead concentrate on the case of a pseudogap DOS, 

piuj) = ^^He{\uj\-D) (12) 

with r > 0. This is the situation of a semimetal with vanishing DOS at the Fermi level. 
Consequently, the tendency toward Kondo screening is reduced, such that no screening 
occurs at small Kondo coupling Jq. As a result, a quantum phase transition between 
phases without and with screening occurs upon increasing Jq [7, 27], as discussed in 
detail below. Importantly, the form of the DOS (12) with r = 1 is relevant for both 
d-wave superconductors and charge-neutral graphene at low energies. 

3.2. Phase diagram 

Despite its simplicity, the pseudogap Kondo model has an extraordinarily rich phase 
diagram, first determined by Gonzalez-Buxton and Ingersent [17] using Wilson's 
numerical renormalization group (NRG) technique. The physics depends not only on 
Jo and the exponent r of the low-energy DOS, but also on the presence or absence of 
particle-hole symmetry. 

Thanks to both numerical [17, 28, 29] and perturbative [30, 31] RG studies, the 
phases and phase transitions of the pseudogap Kondo and Anderson models are by now 
fully understood, and will be summarized below. This discussion, restricted to r > 0, 
mainly follows Ref. [17] and is cast in the language of RG fiows and fixed points. We 
start by listing the fixed points of the problem. For brevity, we will use the acronyms 
of Ref. [17] for both the phases and their stable fixed points. 

Local-moment phase LM. Here the impurity moment is asymptotically decoupled from 
the host and behaves like a free local moment, i.e. it has a residual entropy 5'imp = In 2. 

Symmetric strong- coupling phase SSC. This phase corresponds to Kondo screening in 
the presence of particle-hole symmetry and is the generalization of the metallic Kondo- 
screened phase to finite r. Interestingly, SSC is characterized by a residual entropy 
'S'imp = 2rln2, which means the impurity moment is only partially screened. 
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Figure 3. Schematic RG flow diagrams for the pseudogap Kondo model [17] in 
the plane spanned by the Kondo coupling Jq and the potential scattering Vq, the 
latter measuring particle-hole asymmetry. The flow topology changes qualitatively as 
function of r, as shown in the three panels, with r* = 0.375 ± 0.002 and Tmax = 1/2, 
for details see text. Full dots denote stable fixed points, while open dots are critical 
fixed points. Dashed lines denote separatrices, i.e., phase boundaries. 

Asymmetric strong- coupling phase ASC. In the absence of particle-hole symmetry and 
for r > 0, full screening with Si^p — obtains at ASC, which is maximally particle-hole 
asymmetric. 

Critical points SCR and ACR. The pseudogap Kondo problem has two distinct 

critical fixed points, distinguished by their symmetry under particle-hole transformation 
and denoted "symmetric critical" (SCR) and "asymmetric critical" (ACR), respectively. 
They both exist over a restricted range of r values. 

As deduced from the numerical solution of the pseudogap Kondo model [17], the 
topology of the phase diagram changes qualitatively as the bath exponent r is varied. 
Different phase diagram topologies are observed in three regimes, see Fig. 3. 

a; < r < r* - 0.375 ± 0.002 

(i) For particle-hole symmetry, a critical coupling J^., associated with SCR, separates 
LM from SSC. For initial values J < the flow is directed towards LM, whereas 
for J > Jc the flow is directed towards SSC. 

(ii) For finite particle-hole asymmetry, i.e. Vq 7^ 0, there is a separatrix which separates 
the flow towards LM from the flow towards ASC. 

(iii) Particle-hole asymmetry is irrelevant at LM and SCR while it is relevant at SSC 
where it drives the flow towards ASC. SCR is thus a multicritical fixed point. 

b) r* <r < Frnax = 1/2 

(i) For Vo = there still exists a critical coupling which separates LM from SSC. 

(ii) SCR is now unstable w.r.t. particle-hole asymmetry, and a new asymmetric critical 
fixed point ACR emerges, controlling the transition between LM and ASC. 
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(i) SCR merges with SSC, such that there is no Kondo screening at particle-hole 
symmetry, irrespective of the strength of the Kondo coupling Jq. 

(ii) Screening is possible for finite asymmetry, where ACR continues to control the 
LM-ASC transition. 

(iii) The critical exponents at ACR take trivial values for r > 1, such that r — 1 - the 
case relevant for charge-neutral graphene - acts as an upper-critical dimension. 

Notably, the pseudogap Kondo and Anderson models share identical fixed points 
and quantum phase transitions [17]. This observation can be rationalized within the 
effective field theories described in Sec. 3.4 below. 

3.3. Slave-boson mean- field theory 

A simple and popular approach to the Kondo model in Eq. (9) is the slave-boson mean- 
field theory [2, 32]; a very similar mean-field theory can be applied to the Anderson 
model (5). In this approach, formally justified in a limit where the spin symmetry is 
taken to be SU(A^) with — )■ oo, the Kondo interaction is decoupled by a Hubbard- 
Stratonovich field which is then approximated to be static. This results in a renormalized 
free-particle Hamiltonian which reproduces salient low-temperature properties of a 
Kondo-screened impurity in a metal. 

The slave-boson approach has been applied to the pseudogap Kondo model in 
numerous of papers [16, 27, 33, 34, 35, 36, 37]. It reproduces the existence of a quantum 
phase transition at r > 0, however, the critical properties of this transition are only 
correctly captured for small r. For larger r including the graphene case r = 1, the slave- 
boson method becomes unreliable. It is not sensitive to the subtle effects of particle-hole 
symmetry breaking: It fails to describe the properties near the ACR fixed point, and it 
incorrectly predicts a phase transition for r > 1/2 even in the particlc-holc symmetric 
case. Therefore, quantitative calculations require numerical methods such as NRG. 

3.4- Critical field theories 

The complicated topology of the RG flow, Fig. 3, suggests that different field theories 
are required to describe the critical properties near the SCR and ACR fixed points. 
Such field theories have been worked out in detail in Refs. [30, 31] and provide an 
essentially complete analytical understanding of the pseudogap Anderson and Kondo 
models. Interestingly, none of these field theories is of conventional (i.e. bosonic) 
Landau-Ginzburg- Wilson type; instead all are of genuinely fermionic character and are 
formulated in the degrees of freedom of either the Kondo or the Anderson model. 

In the following we shall summarize the three relevant critical theories. As will 
become clear, only the third will be appropriate to describe the quantum phase transition 
of Kondo impurities in charge- neutral graphene, where r = 1. When specifying flow 
equations from perturbative RG, we will assume a symmetric pseudogap density of 
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states as in Eq. (12). The effect of a high-energy particle-hole asymmetry in the DOS 
can absorbed in the impurity part of the Hamiltonian, e.g., the potential scattering 
term of the Kondo model. This can be rationalized within RG, where integrating out the 
particle-hole asymmetric piece of the bath at high energies yields an effective model with 
particle-hole symmetric bath at low energies and a renormalized impurity Hamiltonian, 
where in particular the particle-hole asymmetry is accumulated. 



3.4- 1- SCR: Kondo model. For small r, an efficient description of the physics at SCR 
is obtained via the Kondo model itself, Eq. (9). A perturbative expansion can be 
performed in Jq and Vq around the LM fixed point where Jq = Vq = 0. As is standard 
practice, we introduce dimensionless couplings j and v, for details see Ref. [31]. Power 
counting reveals that both couplings are marginal for r ~ and irrelevant for r > 0, 
dim[7] = dim[t'] = — r. The one-loop flow equations read 



dlnD '^-^ rfhTD^^"' ('^^ 

where D denotes the running UV cutoff, initially set by the width of the host band. 
Eq. (13) yields a critical fixed point (SCR) at j* — r-l-C(r^), v* = 0, which separates the 
flows towards weak and strong coupling. Controlled calculations near SCR are therefore 
possible in a double expansion in r and j. Potential scattering is irrelevant at SCR and 
consequently does not play a role for leading critical exponents. 

Comparing these properties with the numerically deduced flows in Fig. 3, it is clear 
that this Kondo description of SCR is appropriate for < r < r*. It does, however, 
not capture the physics for r > r* where v becomes a relevant perturbation at the SCR 
fixed point, and it is obviously inappropriate for the graphene case r = 1. 

A perturbative calculation of static critical properties of SCR using the Kondo 
expansion indeed shows excellent agreement with NRG results for small r [31]. Crossover 
functions and dynamical properties have been studied as well, using a combination of 
perturbative RG and Callan-Szymanzik equations [38] , and agreement with NRG results 
has been found where those are available. || 

3.4-2. SCR: Symmetric Anderson model. The Anderson impurity model, originally 
introduced as model for local- moment formation (Sec. 2.2), turns out to provide the 
relevant degrees of freedom to describe pseudogap Kondo criticality for all r > [31]. 
This generically implies that critical fluctuations occur not only in the spin channel, but 
also in the charge channel [40] . 

To discuss the critical behaviour near SCR, we consider a symmetric AIM, Eq. (5) 
with ed — —U/2, a momentun-independent hybridization v, and a particle-hole 

II In the quantum-relaxational finite-temperature regime of a; <i; T numerical studies are notoriously 
difficult. Here, the Callan-Szymanzik results [38] disagree with that from a dynamical large- A'' 
study [39]: the latter indicates that the imaginary part of local Green function (or T matrix) G"'(a;,T) 
multiplied by T'' goes to a non-zero constant, i.e. T' G" (to/T — > 0) = c 7^ 0, while the former suggests 
that it goes to zero. 
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symmetric bath DOS as in Eq. (12). The point ea — U — v — Qis referred to as the free- 
impurity fixed point (FImp), whereas the parameter sets v = and — —U/2 = ±oo 
correspond to doubly degenerate local-moment states in the charge and spin channel, 
respectively. Therefore f = 0, = — oo can be identified with the LM fixed point, while 
V = 0, = cxD is dubbed LM'. 

Notably, the Anderson model is exactly solvable for any v ai U — 0, known as 
resonant-level model. In the particle-hole symmetric case, its low-energy physics can be 
identified with that of the SSC fixed point introduced above:^ its properties correspond 
to a partial screening of the impurity degrees of freedom, with a residual entropy of 



A perturbative expansion is now possible in U around the SSC fixed point. The 
scaling dimension of the renormalized Coulomb interaction u at SSC is found to be 
dim[ti] = — f = — (1 — 2r). The RG flow of u to two-loop order reads [31] 



This flow, together with the trivial flow near LM, LM', and FImp is illustrated in Fig. 4 
[31]. We have included the metallic case, r — 0, Fig. 4a, where the Kondo-screened flxed 
point (SC) is the only stable flxed point. For r > 0, LM is stable as well, while SSC 
becomes unstable for r > Tmax, as can be seen from Eq. (14). Therefore, a critical flxed 
point (SCR) emerges for < r < rmax, Fig. 4b, consistent with Fig. 3. Its properties 
can now be accessed in a double expansion in r and u, and Eq. (14) yields for the 
fixed-point coupling at u*^ = ^^^2in4) ^' ^ perturbative calculation of static critical 
properties again yields excellent agreement with NRG results, here for r < r^ax [31]- 

Owing to particle-hole symmetry, the behaviour at e > is formally identical to 
that at e < 0, Fig. 4, with the latter describing spin-Kondo physics while the former 
corresponds to charge-Kondo physics. Finally, we note that particle-hole symmetry 
breaking can be studied perturbatively for small f, and is found to be relevant at both 
SSC and SCR flxed points, again consistent with Fig. 3. In fact, estimating its scaling 
dimension at SCR as function of r, one flnds that particlc-holc symmetry breaking 
changes from being relevant to being irrelevant once r is reduced below r ~ 0.40 [41] - 
this value can be identifled with r*, Fig. 3. 

3.4-3. ACR: Asymmetric Anderson model. We now turn to the ACR flxed point 
present for r > r*. This is ultimately relevant for understanding the Kondo effect 
in graphene, because there r = 1 at charge neutrality, and t' breaks particle-hole 
symmetry already on the level of the band structure. 

It was realized in Ref. [30] that the critical theory for ACR is that of a level crossing 
of a many-body singlet and a many-body doublet, minimally coupled to conduction 

^ In the metallic case, r = 0, it is known that the resonant-level model is the correct fixed-point theory 
for a Kondo-screened impurity spin. 



imp 



2rln2 [17, 31]. 




(14) 




Figure 4. RG flow diagrams for tlic particle-holc symmetric Anderson model [31], in 
the plane spanned by the level energy e = —u/2 and the hybridization v^. Symbols 
are as in Fig. 3; the flow of u near SSC is in Eq. (14). a.) Metallic case r = 0. SC 
denotes the fixed point of the resonant level model, with conventional Kondo screening. 
LM and LM' correspond to local moments formed in the spin or charge channel, b.) 
< r < 1/2: The critical fixed point SCR (SCR') divides the flow to LM (LM') from 
that to SSC. c.) 1/2 < r < 1: SCR and SCR' merge with SSC as r 1/2", such that 
SSC is now unstable, d.) r > 1: SSC merges with FImp at e = ?; = as r 1~. For 
all r > 1/2, LM and LM' are the only stable phases in the presence of particle-hole 
symmetry. For details see text and Ref. [31]. 



electrons. Using the notation of Rcf. [30, 31], its Hamiltonian can be written as 

^ = 5^ ^Aa^'ka + + 30 [W){s\Ca{^) + h-c] (15) 

where |(t) = I 'I"), I J,) and \s) represent the three allowed impurity states. Sq is the tuning 
parameter ("mass") of the QPT, i.e. the (bare) energy difference between doublet and 
singlet states. The QPT occurs at some sq = Sc, with screening present for Sq > Sc- 
Remarkably, this theory is identical to a maximally particle-hole asymmetric Anderson 
impurity model, Eq. (5), where the doubly occupied state has been projected out, 
U — >■ oo, and (so^S'o) in Eq. (15) have been identified with (e^, t;) in Eq. (5). 

In this model, the point £o = S'o = is dubbed valence-fluctuation fixed point 
(VFl). As above, Qq = 0, £o = "OO corresponds to LM, while = = oo describes 
a fully screened and particle-hole asymmetric singlet state, to be identified with ASC. 

A perturbative expansion is now possible in qq around VFl. Power counting yields 
the scaling dimension of the renormalized hybridization dim[(7] = f = The one-loop 
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a) b) 

Figure 5. RG flow for the asymmetric Anderson model [30, 31] in the e-g"^ plane, 
obtained from Eq. (16) a.) r < 1: The critical fixed point ACR separates the flow 
towards LM from that to ASC. b.) r > 1: ASC merges with VFl as r — >■ 1~, which 
describes a level crossing with perturbative corrections. The behaviour near r = 1 
is similar to that of the Landau-Ginzburg model near d = 4, with VFl and ASC 
corresponding to the Gaussian and Wilson-Fisher fixed points, respectively. 



flow equations for g and the renormalized mass s read 



dg _ 3 



3 



^=-.-/ + 3/.. (16) 

results to two-loop order can be found in Ref. [31]. The RG flow is shown in Fig. 5 - 
this flow has strong similarity to that of the standard Landau-Ginzburg model. The 
fact that g is relevant for r < 1 and irrelevant for r > 1 allows us to identify r — 1 
as an upper-critical dimension of the pseudogap Kondo problem, akin to (i = 4 in the 
Landau-Ginzburg theory. For r < 1 a non-trivial fixed point (ACR) emerges at g*^ = |f 
and e* = —'^r, Fig. 5a, similar to the celebrated Wilson- Fisher fixed point. Critical 
properties, evaluated in a double expansion in f and g, again agree well with NRG 
results [31]. In contrast, for r > 1 in Fig. 5b, we have "Gaussian" behaviour controlled 
by the VFl fixed point, which here corresponds to a simple level crossing with corrections 
captured by plain perturbation theory in go. In the case r — 1, relevant to charge-neutral 
graphene, this perturbation theory is logarithmically divergent at criticality and needs 
to be resummed, as is standard at the upper critical dimension. 

The structure of the critical theory (15) implies that spin and charge fiuctuations are 
strongly coupled, i.e. suitably defined observables in the charge sector become critical 
at the Kondo quantum phase transition controlled by ACR [30, 40]. 

3.5. Finite carrier concentration 

For graphene away from charge neutrality, yU 7^ 0, the DOS at the Fermi level is finite, 
and consequently a magnetic impurity described by the Kondo model (9) will be screened 
in the low-temperature limit for any value of the Kondo coupling Jq [42, 44, 43]. To be 



Kondo Impurities in Graphene 



16 




(17) 



In the limit of small Jq the corresponding Kondo temperature will be exponentially 
small according to InTx oc but in general Tk needs to be calculated numerically, 

as the simple equation (11) is no longer applicable due to the strong energy dependence 
of the DOS. 

On general grounds, one expects that the presence of the quantum phase transition 
at Jo = Jc, = influences the behaviour at finite // as well. Fig. 6a. In this quantum 
critical regime, heuristic scaling arguments suggest Tk = k,\ii\ with a universal constant 
K depending on r only. This problem can be tackled by generalizing the RG equations 
(16) obtained for ACR to a finite chemical potential [43]: 



with Fi^2{y) — |l + yr±2|l— yl''. The last term in de /din D describes the level shift due 
to the real part of the bath Green's function. A detailed discussion of these equations 
has been given in Ref. [43], showing that the asymmetric nature of the critical theory 
induces a strong asymmetry between the two signs of /i in Tk(/x). In fact, negative 
yU drives the near-critical system directly into the screened phase, whereas positive n 
first induces a crossover to a spin- 1/2 moment which is subsequently screened via a 
conventional Kondo effect. For r < 1, where the ACR fixed point is interacting, the 
scaling prediction holds with Tk = ft;±|/i| for ^ 0. In contrast, at the upper critical 
dimension, r — 1, Tk = continues to hold for /i < 0, while for > logarithmic 

corrections and Kondo logarithms conspire such that Tk oc |//|^ where x ~ 2.6 is a 
universal exponent. 

This quantum critical particle-hole asymmetry of Tk(/x) pertains to the off-critical 
situation as well. This is nicely seen in the numerical results. Fig. 6b, obtained for a 
realistic graphene DOS: There is not only asymmetric behaviour for Jq < Jc, but also 
the minimum of Tk(/x) for Jq > Jc is not found at = but somewhat away from it. 

3.6. Practical consequences for Kondo screening in graphene 

Based on the analysis presented so far, we can now specify theoretical predictions for the 
screening of Kondo impurities in graphene. The half-metallic energy dependence of the 
graphene DOS implies that Kondo screening tends to be weaker than in conventional 
metals, implying smaller Tk, and that the single-parameter scaling known from the 
metallic Kondo problem [2] will in most cases not apply. Such unconventional multi- 
scale crossovers have been studied explicitly [43, 44]. 

The absence of single-parameter scaling may in fact complicate the unambiguous 
identification of the Kondo effect - note that typically the scahng of the susceptibility 




(18) 
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Figure 6. a) Phase diagram for the pseudogap Kondo problem for r = 1 in tlie 
presence of finite chemical potential /i. A quantum phase transition between LM and 
ASC exist in the charge-neutral case, ^ = 0, whereas the local moment undergoes 
Kondo screening for all /i ^ 0. The shaded region is influenced by quantum criticality, 
with T'K(/i) following a power law, for details see text and Ref. [43]. b) NRG 
results [43] for the Kondo temperature Tk sls function of /x for different values of 
the Kondo coupling Jq, calculated for a DOS appropriate for Co on graphene [12] 
where X. « 4.3 eV. 



or the resistivity correction has been used to argue in favor of Kondo screening. We 
beheve that, instead, the doping dependence of Tk, the latter extracted, e.g., from the 
peak width of a tunnehng spectrum, can be used as a key indicator of Kondo screening: 
Tk will strongly vary with doping, with a minimum near charge neutrality. 

A concrete prediction of Tk requires a microscopic modelling for the specific type 
of impurity, which would yield the hybridization and interaction terms of an Anderson 
impurity model. As discussed in Sec. 2, at present there are considerable uncertainties 
in these parameters for all relevant impurities. Nevertheless, one may use the ab-initio 
results for the hybridization function of a Co adatom with S* = 1/2 in an h position 
[12] (which should apply to NiH as well) to numerically calculate Tk within an effective 
Kondo model. Results obtained using NRG have been presented in Ref. [43] and are 
reproduced in Fig. 6b. The strong fi dependence and the pronounced electron-hole 
asymmetry of Ti<i{fi) are apparent. At Jq = Jc the linear and power-law behaviours of 
TKifJ'), advertised in Sec. 3.5, are nicely visible. Comments on numbers are in order: 
(i) Due to the uncertainty in the Coulomb interaction U, the effective Kondo coupling 
Jo is not known to a good accuracy. Therefore, any prediction of Tk, in particular near 
charge neutrality, comes with excessively large error bars, (ii) In the model leading 
to Fig. 6b, we have assumed SU(2) symmetry. However, according to Ref. [12], the Co 
impurity has an approximate SU(4) symmetry which is broken down to SU(2) on a scale 
of 60 meV. Therefore, the high-energy flow of the Kondo coupling will differ for the two 
models, such that the critical coupling for Co in this S* = 1/2 state is predicted to 
be 2.2 eV, approximately matching the estimate in Ref. [12]. 
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3.7. STM and quasiparticle interference 

Conduction-electron scattering off impurities can be probed using STM: Conductance 
spectra can be recorded at/near the impurity site and as function of the distance to the 
impurity. In particular, spatial variations in the local DOS, p(f, uu), can be interpreted in 
terms of impurity-induced energy-dependent Friedel oscillations, so-called quasiparticle 
interference (QPI). Analyzing QPI spectra using models of elastic scattering allows to 
extract information on both the host band structure and the nature of the impurity. 

For graphene, the initial experiments of Manoharan [45] have triggered a number 
of theoretical studies of local spectra [36, 44, 46, 47] and the expected QPI signal 
[12, 48, 49]. A striking feature of the Kondo effect in charge-neutral graphene is that 
impurity spectral density is not peaked at the Fermi level, but away from it, with 
vanishing spectral weight at Ep [44, 50]. In the doped case, single-parameter scahng is 
again violated for a large range of parameters [43, 44]. We note, however, that a detailed 
comparison between theory and experiment is lacking to date. 

4. Experiments and open issues 

Despite numerous attempts to create and study Kondo impurities during several years 
of graphene research, the amount of experimental data is still somewhat limited. Both 
magnetic adatoms and vacancies have been considered, but clear-cut observations which 
unambiguously verify available theories have not been reported to our knowledge. We 
start by summarizing the most prominent experiments, and then discuss issues which 
might contribute to complicate the interpretation of the data. 

4.1. Experiments 

4.1.1. Adatoms. Isolated Co adatoms placed on top of a graphene sheet have been 
studied using STM in Refs. [45, 51]. Ref. [45] employed a conducting SiC substrate 
which shifts the chemical potential to roughly fi = 0.25 cV. Spectral signatures of Kondo 
screening were observed, with Tk ~ 15 K, including the expected splitting of the Kondo 
peak upon apphcation of a magnetic field. Surprisingly, Tk was almost identical for Co 
atoms placed on the t and h positions of a carbon hexagon. Fig. 2a. For Co in the h 
position, the energy dependence of the small-bias conductance was interpreted in terms 
of a two-channel Kondo (2CK) effect, whereas single-channel Kondo (ICK) behaviour 
was found for Co in a t position. 

In contrast, Ref. [51] worked with an insulating Si02 substrate which allowed 
for gate tuning of the chemical potential. Here, clear-cut Kondo signatures were not 
observed, instead the spectral features mainly reflected charging effects and vibrational 
excitations. 

4.1.2. Defects. In a second group of experiments, point defects within the graphene 
sheet were created either by irradiation [52, 53], with estimated defect densities ranging 
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Figure 7. Kondo temperature Tk for vacancy moments in irradiated grapiiene for two 
different samples as function of the gate voltage which controls the carrier concentration 
[52]; V^^min corresponds to charge- neutral case, i.e. the chemical potential located at 
the Dirac point. The Tk values have been extracted from the temperature dependence 
of the electrical resistivity (taken from Ref. [52]). 

from 10~^ [52] to 0.1 [53] per C atom, or by depositing fluorine adatoms [53, 54]. In the 
case of irradation, it is believed that the main defects are carbon vacancies. 

Ref. [52] studied magnetotransport through an irradiated graphene sample placed 
on a Si substrate and found a resistivity increase at low temperatures combined with 
negative magnetoresistance, consistent with Kondo screening. Using extensive fits of 
p(T) to the standard theory of the metallic Kondo effect, Tk was found to vary between 
30 and 90 K for estimated chemical potentials within [-0.3eV,0.3eV], see Fig. 7. Such a 
Tk variation appears rather small compared to that expected within a pseudogap Kondo 
model, see Sec. 3.6 and Fig. 6. Also, it is surprising that the transport data could be 
fitted to theoretical results for metallic Kondo screening even near neutrality where 
standard single-parameter scaling is not expected due to strong energy dependence of 
the host DOS. Despite interesting proposals and ideas [55, 56], we feel that a convincing 
theoretical explanation for the results of Ref. [52] is missing to date. 

Ref. [53] aimed at quantifying defect-induced magnetism in graphene laminates, 
utilizing magnetization measurements away from the dilute limit. For both fiuorine 
adatoms and vacancies, paramagnetic behaviour of spin-1/2 moments was detected. The 
measured magnetic moment per defect was between 0.1 and 0.4 /x^ in the vacancy case, 
possibly consistent with one spin-1/2 per vacancy in the dilute limit (due to uncertainties 
in estimating the vacancy density). For fiuorinated graphene the magnetic moment 
per adatom was only 10"'^ /i^ - this extremely small value was ascribed to adatom 

+ The Kondo interpretation of the transport data of Ref. [52] has been questioned in Ref. [57], where 
it was instead proposed that electron-electron interactions in the presence of disorder are responsible 
for the logarithmic resistivity increase at low temperature. 
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clustering, such that only larger clusters contribute one spin-|. In these experiments, 
neither magnetic order nor signatures of Kondo screening were detected down to 2 K. 

Finally, Ref. [54] investigated weak-localization physics in gated graphene with 
fluorine adatoms with dilute concentrations of order 10~^. The results of the transport 
measurements were interpreted in terms of suppressed weak localization due to spin-flip 
scattering from fluorine-generated moments, but the Kondo temperature was estimated 
to be as small as 0.01 K for carrier densities of 0.6 x 10^^/cm^. Such a low Tk would imply 
a very small magnetic coupling, J ~ 5 meV, between the fluorine-induced moments and 
the Dirac electrons of graphene, possibly consistent with the absence of both magnetic 
order and screening as measured in Ref. [53]. 

4-2. One-channel vs. two-channel Kondo screening 

The two-channel Kondo (2CK) effect emerges if a magnetic impurity is coupled 
symmetrically to two equivalent screening channels of conduction electrons, such that a 
standard Kondo singlet is unstable [58]. Instead, the low-T behaviour is then governed 
by a non-trivial intermediate-coupling flxed point with non-Fermi-liquid properties. 

2CK physics being relevant for graphene has been proposed theoretically in 
Refs. [37, 42, 59], based on the idea that the electrons in the two valleys, i.e., near K and 
K', could form independent screening channels. The resulting pseudogap 2CK model 
displays an interesting interplay of pseudogap Kondo physics and the non-Fermi liquid 
behaviour of the 2CK effect and has been studied in Refs. [17, 60]. However, an analysis 
of possible microscopic models for graphene impurities suggests that unavoidable inter- 
valley scattering will invariably couple the two screening channels such that single- 
channel Kondo screening prevails at least at low energies and temperatures [59], 
i.e. below a crossover scale Tick- For well localized magnetic moments, inter- valley 
scattering is strong resulting in Tick ~ Tk, such that there is unlikely to be an 
intermediate regime of 2CK screening. The interpretation of STM data in terms of 
a 2CK effect in Ref. [45] is therefore puzzling. 

4-3. Influence of electron-hole puddles 

A serious complication for quantitative comparisons between theory and experiment, 
in particular for graphene sheets near charge neutrality, is the presence of electronic 
inhomogeneities, known as electron-hole puddles [61, 62]. These puddles can be 
understood as spatial variations of the local carrier concentration or, equivalently, of 
the local Dirac-point energy, with the characteristic length scale (or puddle size) of 
10-20 nm. While the origin of these puddles is under debate [62, 63, 64], a plausible 
explanation is the influence of charged defects in the substrate. 

In the context of Kondo impurities, these puddles imply a spatial variation 
of the local DOS, which itself determines the local Kondo temperature. Thus, a 
distribution of Kondo temperatures will be present in a sample with multiple impurities. 
Importantly, different macroscopic observables will be dominated by different parts of 
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this distribution: Whereas the magnetic susceptibihty will be dominated by weakly 
screened (i.e. 1ow-Tk) moments, the electric resistivity receives mainly contributions 
from strongly screened (i.e. high-Tx) moments. For metallic systems with a broad 
distribution of Kondo temperatures, non- Fermi-liquid behaviour may arise, as has been 
discussed in the framework of so-called Kondo-disorder models [65, 66]. 

For graphene, an interesting question is whether Kondo disorder could explain 
the weak gate- voltage dependence of the transport-TK in the experiment of Ref. [52]. 
However, we beheve that this is unhkely to be the case: Weak disorder is insufficient to 
significantly modify the Tx{iJ>) dependence, whereas strong disorder would yield a broad 
Tk distribution which appears incompatible with the fact that the transport data [52] 
could be well described by the universal metallic Kondo behaviour with a single Tk- 

4.4- Multiple impurities: RKKY interaction and magnetic order 

In samples with a non- vanishing concentration of magnetic impurity moments, the issue 
of impurity- induced order becomes relevant, with the coupling between the moments 
mediated by conduction electrons via the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
interaction. For magnetic moments in metals, magnetic order sets in only at rather 
high concentrations (typically 20% or more), because otherwise Kondo screening is likely 
to prevail. This can be different in semiconductors, and indeed semiconductors doped 
with magnetic ions display ferromagnetism with Curie temperatures above 100 K [67] , 
possibly important for spintronics applications. 

For graphene, the RKKY interaction has been discussed theoretically [68, 69, 
70]. For charge- neutral graphene, the oscillations with distance r typical of RKKY 
interactions are absent, such that the interaction is strictly ferromagnetic for moments on 
the same sublattice and antiferromagnetic for moments on different sublattices, in both 
cases falling off as Thus, an inhomogeneous but unfrustrated antiferromagnetic 

ground state can be expected for randomly placed moments. Such order is predicted to 
persist at finite temperature even in this 2d situation due to the long-range nature of 
the RKKY interaction [71] . Departing from charge neutrality, a crossover to standard 
metallic behaviour is expected. 

The competition of RKKY-induced order with Kondo screening has not been 
studied in detail, but it is plausible that most considerations for magnetic moments 
in metals apply: Away from charge neutrality, one expects a quantum phase transition 
between a Fermi liquid with screened moments and an ordered state which is either 
antiferromagnetic or spin-glass-like. As with other quasi-2d metals, the nature of this 
phase transition is an open problem [72, 73, 74]. In the neutral case, the transition is 
between a semimetal and a magnetic insulator - such a transition has been investigated 
in the absence of quenched disorder [75], but the effect of randomness due to moment 
disorder has not been studied. Experimentally, impurity-induced magnetic order in 
graphene has not been observed to our knowledge (see e.g. Ref. [53] for an attempt). 
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4-5. Summary 

Studying magnetic moments in graphene holds the prospect of observing exciting 
phenomena, such as Kondo physics beyond one-parameter scahng, single-impurity 
quantum criticality and associated local non- Fermi-liquid behaviour, magnetic order 
from dilute moments, and lattice quantum criticality in strictly two dimensions. While 
for most of these, a theoretical framework is available, experimental data is scarce and 
more experimental activities are clearly called for. 

5. Outlook: Beyond graphene 

The Kondo effect in graphene is related to a variety of other quantum impurity problems. 
This section will highlight the most important connections. 

5.1. Kondo impurities in unconventional superconductors 

The pseudogap Kondo problem was first discussed [27] in the context of magnetic 
impurities in unconventional superconductors. For BCS states with nodes in the 
superconducting gap function, the density of states of Bogoliubov quasiparticles vanishes 
algebraically as u ^ 0. Moreover, a locally coupled impurity, Vk = V in Eq. (5), 
is not directly influenced by pairing, because the local anomalous Green's function 
vanishes in an unconventional superconductor [76]. Then, the Kondo problem in such 
a superconductor is equivalent to the pseudogap Kondo problem described in Sec. 3, 
with a \u\^ density of states at low energies and r — 1 {r — 2) for d-wave (p-wave) 
superconductors. Importantly, the Fermi energy is always pinned to charge neutrality, 
i.e. the crossovers described in Sec. 3.5 cannot be accessed. 

A clear-cut experimental observation of Kondo screening in unconventional 
superconductors ideally requires Tk ^T^. with being the superconducting transition 
temperature, together with the ability to tune either the host gap or the Kondo coupling. 
While the latter can be efficiently varied in nanostructurcs with quantum dots, we are 
not aware of corresponding experimental results, likely because of the lack of suitable 
unconventional superconductors as lead materials. In bulk superconductors with 
impurities, signatures of Kondo screening have been detected in NMR measurements 
on Zn-doped high-7; cuprates of the YBa2Cu302/ family [77, 78]. Similarly, STM 
experiments on Bi2Sr2CaCu208+5 have detected large low-energy conductance peaks 
near Zn impurities [79], which subsequently have been interpreted in terms of a Kondo 
resonance [34, 50].* It should be noted, however, that other interpretations for the STM 
data have been put forward as well [79, 80, 81], and a concise picture for impurity effects 
in cuprates has not yet emerged. 

* For cuprates, it has been assumed that the non-magnetic Zn ion induces a magnetic moment by a 
mechanism of dimer breaking, as happens in spin-gapped magnets such as spin ladders. 
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5.2. Kondo impurities on the surface of topological insulators 

The 2d surface states of 3d topological insulators (TI) admit a low-energy description 
in terms of a Dirac equation. The resulting electronic properties are therefore similar 
to that of graphene, with a few important differences: (i) there is a single Dirac cone 
(or, more generally, an odd number) per surface, and (ii) the role of the pseudospin 
(or sublattice) in the graphene case is taken by the physical spin, such that TI surface 
states display spin-momentum locking, and there is no additional spin degeneracy. 

The physics of Kondo impurities in this setting has been analyzed theoretically 
in a number of papers recently [82, 83, 84]. The main conclusion is that, despite the 
non-trivial topological structure of the TI surface states, the corresponding local Kondo 
problem for a spin S — 1/2 impurity can be mapped onto the standard pseudogap Kondo 
model of Sec. 3. Consequently, strong deviations from conventional, i.e. metallic, Kondo 
screening are expected once the Fermi level is tuned close to the Dirac-point energy of the 
surface states. Experimentally, magnetic moments on TI surfaces have been investigated 
in a few papers [85, 86, 87], but not in the dilute regime with focus on Kondo screening. 

5.3. Kondo physics in spin liquids 

A somewhat different, but still overlapping, topic is the physics of magnetic impurities 
embedded in quantum magnets. The closest relation to the graphene Kondo problem 
is found for host magnets without scmiclassical long-range order, but still a gapless 
spectrum of excitations. Two interesting cases will be discussed in the following. 

In quantum-critical magnets, located near a zero-temperature transition between an 
antiferromagnet and a paramagnet, the elementary host excitations are spin-1 critical 
magnons which interact with the impurity spin via a Yukawa-type (i.e. three-point) 
coupling. As a result, true Kondo singlet formation is not possible. In dimensions 
1 < d < 3, RG studies have predicted a partial screening of the spin, described by a 
non-trivial intermediate-coupling fixed point [88, 89]. In d — 2 this prediction has been 
verified numerically [90]. 

Magnetic impurities in gapless spin-liquid phases are more akin to Kondo problems: 
Here the host excitations arc typically spin-1/2 spinous coupled to a U(l) and Z2 gauge 
field. The spinous couple to the impurity spin with a Kondo-like (i.e. four-point) 
interaction. The physics depends on the nature of the spinous, and a few cases have 
been discussed in the literature. Linearly dispersing bosonic spinous yield a rich phase 
diagram, with a variety of possible T — phases, including the possibility of full Kondo 
screening, and quantum phase transitions [38]. In contrast, fermionic spinous lead to 
physics similar to standard Kondo expectations: In the presence of a spinon Fermi 
surface, the impurity spin gets always screened at low T [91]. In the case of 2d Dirac 
spinous of an algebraic spin liquid, a quantum phase transition not unlike that of the 
pseudogap Kondo problem. Sec. 3, emerges [33, 92, 93]. However, it should be noted 
that the infiuence of gauge fields beyond perturbation theory has been neglected in the 
published treatments. 
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